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We present results of an ongoing study of two flavour QCD with Wilson fermions at finite temperature. We 
have used tree level Symanzik improvement in both the gauge and fermion part of the action. The phase diagram 
was previously determined on an 8^ x 4 lattice and the existence of Aoki's phase demonstrated. On our current 
lattice of 12^ x 24 x 4 we have extended the set of observables and studied chiral and thermal properties at light 
quark masses. 



1. Introduction 



2. Results 



The study of the finite temperature phase 
diagram of 2 flavour QCD with Wilson quarks has 
revealed a rather intricate picture. As was pointed 
out already some time ago [Q, the pion can 
become massless despite the fact that the Wilson 
term breaks chiral symmetry, because of the 
existence of a second order phase transition into 
a parity and flavour symmetry violation phase. 
At finite temperature this phase is expected to 
pinch into a cusp that meets with the m — line 
coming in from weak coupling The thermal 
line of the deconfinement phase transition crosses 
the rUq = line and in it was argued that it 
should bounce back towards weaker coupling due 
to a symmetry under the change of sign of the 
mass term in the continuum theory . The relation 
of this crossing point to the tip of the cusp of the 
Aoki phase is somewhat subtle. First of all the 
thermal line should not run into the Aoki phase, 
since this would mean, that one has a massless 
pion in the deconfined phase. The thermal line 
can now either touch the cusp of the Aoki phase or 
not. In the latter case one would have a transition 
line (probably first order) from a confined phase 
with niq > to a confined phase with niq < 0, 
without the pion actually becoming massless M. 



Utilizing a tree level 0{a^) improved action for 
the gauge fields and the tree level clover action 
for the fermion fields, we have, in a previous 
study on an 8"^ x 4 lattice, mapped out the phase 
diagram using the pion norm and the Polyakov 
loop as observables. In the present study we have 
increased the lattice size to 12^ x 24 to check finite 
size effects and to measure the pion screening 
mass and quark mass accurately. From these we 
can also determine the properly subtracted order 
parameter of chiral symmetry breaking for Wilson 
fermions. 

2.1. Quark mass 

The current quark mass is defined via the axial 
Ward identity ||] : 



2m„ 



(0\PM 



(1) 



*Work in part supported by NATO grant no. CRG940451 
and by the European Community TMR Programs TRACS 
and ERBFMRX-CT97-0122 



Our results are shown in figure 1. For f3 = 2.8 
we find some curvature for the quark mass as 
a function of 1/k, though a linear fit produces 
a reasonable x^- For (3 = 3.1 we have also 
explored a region of hopping parameters where 
the quark mass becomes negative. There we find 
a rather peculiar behaviour. The quark mass 
decreases as one lowers k towards Kc (see section 
2.3), which is in contrast to simulations with 
unimproved Wilson fermions that have shown 




the same qualitative behaviour for positive and 
negative quark masses |2] . 

2.2. Pion screening mass 

We have also accurately measured the pion 
screening mass. Our results are depicted in figure 
2. For /? = 2.8 the decrease of the pion mass is 
consistent with a linear behaviour oc 1/k — 
1/kc down to small values of . This also applies 
for (3 = 3.1 for n values that correspond to 
positive quark masses. For negative quark masses 
the behaviour is quite different and inconsistent 
with a linear behaviour for the points measured. 
We will argue below that for /3 = 3.1 the pion 
mass does not go to zero as the quark mass goes 
to zero, because one crosses the finite temperature 
transition line before the quark mass becomes 
zero. 

2.3. Locating the critical line 

To locate the critical hopping parameter Kc we 
have performed the following fits to our data. 
We have fitted linearly and niq quadratically 
as a function of Since the pion norm near 
Kc is expected to behave like 11 « l/m^^ we 
have also extracted a value for Kc from this 
observable. Since the data show quite some 
nonlinear behaviour as a function of we have 
used a quadratic fit ansatz. As the results from 
the different fits turned out to be consistent, we 
took confidence in using our data for the pion 



norm from our simulation at (3 = 3.0 on the 
smaller lattice to get an estimate of Kc for (3 = 3.0. 
Our results are summarized in the following table: 
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(3 =3.0 


/?=3.1 




0.1860(2) 




0.1783(2) 


TUq 


0.1853(3) 




0.1770(3) 


n 


0.1859(3) 


0.1823(10) 


0.1800(5) 



2.4. Chiral order parameter 

Because of the explicit breaking of chiral 
symmetry by the Wilson action, one has to 
define a properly subtracted order parameter to 
obtain the correct continuum limit. Using axial 
Ward identities the order parameter is defined as 
follows d): 

(V;V')s«6 = 2mg-Z-^(7r(a;)^(0)} (2) 

Here Z is a renormalization factor for which we 
take its tree level value Z = (2k) ^. The sum 
over the pion correlation function is just the pion 
norm. For (3 = 2.8 the data extrapolate to a finite 
intercept at ruq = 0. For (3 = 3.1 the data show 
more curvature and we expect {ipip)sub to go to 
zero, though a definitive statement has to await 
further data at smaller quark masses. 

2.5. Polyakov loop 

Figure 4 shows our results for the Polyakov 
loop as a function of k including data from both 




Figure 3. Subtracted chiral condensate as a Figure 4. Polyakov loop as a function of k for 
function of 2m, for /3 = 2.8 and (3 = 3.1. /? = 2.8, /3 = 3.0 and f3 = 3.1 



lattices. The vertical lines indicate the position 
of the extrapolated Kc, which decreases as f3 
increases. For f3 — 2.8 all our data for the pion 
mass lie in the confined region. We therefore 
conclude that the pion mass vanishes as the quark 
mass goes to zero, i.e. for /3 = 2.8 we hit the Aoki 
phase as we increase k. For (3 = 3.0 and 3.1 this 
is no longer so clear. The Polyakov loop already 
shows a high temperature behaviour where the 
extrapolated pion mass would be small. For /3 = 
3.0 the situation is less prominent and one could 
still argue that the pion becomes massless, but 
since the Polyakov loop is in the high temperature 
phase immediately after one crosses Hc, the finite 
temperature line and the line of vanishing quark 
mass come very close together. For f3 = 3.1 
it becomes evident, that one crosses the finite 
temperature transition line before the line of 
vanishing quark mass. The pion will therefore 
not become massless as the quark mass vanishes. 
Because the Polyakov loop continues to rise past 
the niq = line, one can exclude that the finite 
temperature line bounces back towards weaker 
coupling immediately. 

3. Summary and Conclusions 

From the vanishing of the pion mass at zero 
quark mass we conclude that at /3 = 2.8 there 
exists a parity and flavour symmetry breaking 



phase. At /3 — 3.0 the finite temperature 
transition line and the line of vanishing quark 
mass come at least very close and for P — 3.1 
the Aoki phase ceases to exist and one crosses 
the finite temperature transition line before the 
line of vanishing quark mass and therefore the 
pion does not become massless. To make these 
claims even stronger it would be desirable to 
lower the quark mass further for /? = 3.1 and 
to see the pion mass rising again. Once the phase 
diagram is known one can start to investigate the 
thermodynamic behaviour at small pion masses. 
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